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1. Introduction. Browsing through an integral table on a dull Sunday 
afternoon some time ago, I came across four divergent trigonometric inte- 
grals. (See (0) and (g) below.) I was intrigued as to how these divergent 
integrals ended up in a respectable table. Tracing their history, it turned out 
they were originally "evaluated" when some convergent integrals, (|5[) and 
(H), were differentiated under the integral sign with respect to a parameter, 
formally yielding (|1|) and (fj). We will give a simple proof that these inte- 
grals diverge, look at their history in print and then make some final remarks 
about necessary and sufficient conditions for differentiating under the inte- 
gral sign. We have no motive in defaming either the (shockingly famous) 
mathematician who made the original error, or the editors of the otherwise 
fine tables in which the integrals appear. We all make mistakes and we're not 
out to point the finger at anyone. (In this regard see the last two exercises 
of Chapter 2 in JHJ.) We will also see that Maple and Mathematica have 
considerable difficulties with these integrals. 

2. Four divergent integrals. Here they are. Throughout, a and b are 
positive real numbers. Purported values appear on the right: 




x=0 

i_ rn b_ nr 

2a \ J T 2a \ 2a 
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The two Fresnel integrals are 

^/ n 1 r dt i n, s i r . dt 

C(x) = , / cosi^= and Six) = , / sint— . (3) 



Note that in the literature, the same symbols C and S denote several different 
definitions of the Fresnel integrals. 
Let's prove these devils diverge. 

Proposition: The integrals in ([[]) and (0) diverge. 

Proof: Consider A := f^° oo xe l( - x2+x ^ dx (which, unfortunately, does not 
exist). Since the integrand is continuous, this integral exists if and only if 
the limits 

lim / xe l( - x2+x ^ dx lim / xe i( - x2+x ^ dx 



T >00 /q T >0O J rp 



both exist. Let T 1; T 2 > 0. Integrate by parts and complete the square: 

T 2 i/4 T2+1/2 

/ 2 xe< x2+ ^ dx = - \e l ^ 2+T ^ - e^ T2 - T A - — / e lx " dx. (4) 

-Ti+l/2 

Let's look at the convergence of / := f °° e" 2 dx. Use the substitution x 2 = 
i. Then / = e^dt/y/i = \/2n(C(oo) + iS'(oo)). This can be seen to 
converge by applying Dirichlet's Test over (l,oo). (For Dirichlet's Test see, 
for example, |p6fl , pp 261.) More properly, we should start with the version 
of /, show it converges and then transform back to the 'x' version. But let's 
keep it cool. In fact, many ways have been found to evaluate I. One method 
is to use contour integration. (Rotate the integral J^° oo e~ x dx by 7r/4 in the 
complex plane. See p9| , pp 184.) A second method is to use the gamma 
function |26j, pp 272. The result is I = e in/4 ^. Now, as T u T 2 -> oo 
(independently of each other) the final integral in (|J) becomes / but the 
bracketed term fails to have a limit. Hence, the integral A diverges. 

To get the integrals in ([I]) and (|2]) we do the following. Suppose B := 
xe 1<yX ~ x ' dx converged. The transformation x i— > — x gives B = —A 
so B diverges. The transformations x i— > x/y/a± (y/a — b)/(2a) now show 
that C := xe l ( ax ±bx ^ dx diverges for all positive a and b. Finally, if the 
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integrals in ([!]) and (0) converged then we could form the four convergent 
linear combinations 



x [cos(ax 2 ) cos(6x) =F sin(ax 2 ) sin(fcr)] dx 

x=0 



oo 



sin(ax 2 ) cos(6x) ± cos(aa; 2 ) sin(for)] dx 



x=0 

But, using the addition formulas for the sine and cosine functions and then 
Euler's formula yields C. Hence, the integrals in (|I|) and (fj) diverge. ■ 
To see the manner in which the integrals diverge, let 

A T := / xe i{x2+x) dx 



T 

e-^ 4 r T+1/2 2 

= e lT sinT-- / e™ dx. 

2 J-T+l/2 

As T — > oo, the integral term in the line above has limit y / 7re^ 7r_1 ' ) / 4 /2, 
whereas the term e tT sinT oscillates rapidly with unit magnitude. Note 
that this also shows that our integrals do not even exist as Cauchy principal 
values. 

3. History of the divergent integrals. Now we'll look at the history 
of (0) and (HP in print. The thickest book of integrals (3500 pages in five 
volumes) is that of PrudnikovQ, Brychkov and Marichev, |25| . Our integrals 
appear in Volume I, 2.5.22. They also appear in the original Russian edition 



24j . Sources are not referenced in this work. It is interesting that they do not 



appear in the earlier book || by Ditkin and Prudnikov. Other major tables 



they are absent from are [0, |]17| , |22j and |23|]. As the tables by Erdelyi 
and Oberhettinger are quite comprehensive, one suspects it was noticed that 
these integrals diverged and they were purposely omitted. However, they 
are contained in the Gradshteyn and Ryzhik tome, |15| (3.851). They are 
not in the first few Russian editions but the 1963 edition Jl3| was enlarged 
enough to include these divergent integrals. All subsequent Russian editions 
and English translations beginning 1965 |14| contain the integrals in (0) and 



^Recently deceased. A heartfelt obituary by Marichev appears in Ell. See also plj 
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(fj). Now, Gradshteyn and Ryzhik do give references. They say (in a garbled 
citation) that our integrals come from tables by Bierens de Haan. 

David Bierens de Haan (1822-1895) was a Dutch mathematician noted 
for compiling tables of integrals, for actuarial work, for writing various essays 
in the history of science and mathematics, for producing an encyclopaedic bi- 
ography of Dutch scientists and for being an early editor of the works of 
Christian Huygens. (A mammoth task. It took until 1950 when the 22 nd 
volume was finally published.) A complete list of de Haan's publications is 
given in f[9| . There have been several papers on his life and work. See |27] 



for references, photos and an interesting reproduction of the 1935 title page 
from a Japanese edition of his integral table. His 1858 Tables d'integrales 
definies |2| was the first really substantial table of integrals. It was enlarged 
and corrected in an 1867 edition Q . For nearly a century these were the pre- 
eminent integral tables. The 1867 edition was still being reprinted in 1957 



||, three years after the publication of the Bateman Manuscript tables |T0 
The integrals in ([!]) appear in the 1858 table f^[, formulas 193.17 and 193.18, 
an 1862 companion volume that details the techniques used to compute in- 
tegrals in the tables [Q, pp 443, and in the 1867 table 0], formulas 150.4, 
150.7. 

Now, Bierens de Haan lists Cauchy as his source for (|l|). An examination 
of Cauchy's works (27 volumes!) shows these integrals appear twice || (1815) 
and || (see also 0) (1825). In both instances, Cauchy correctly obtains the 
convergent integrals 



oo 

sinfaa; 2 ) , , 

2 \ > cos(taj ax 



cos ax 




2 = 



cosi i) Tsin (y 



(5) 



with a = 1. As is clearly stated in the above references, he then proceeds 
to differentiate under the integral sign with respect to b. Very bad! The 
functions defined by the integrals in (Q) are certainly differentiate since the 
right side of (|5]) is differentiable. But, differentiating under the integral sign 
leads to our divergent integrals ([!]). It is not some unknown schmo but 
Cauchy, the "Father of rigour" , who commits an error that has been copied 
for 185 years. (The appropriateness of this epithet is contested. One in 
favour is [12]. One against is |16|.) Bierens de Haan repeats this argument 



in i,pp 443. 
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When the two integrals 

sm(bx) dx 



sin (ax 2 ) 

■2\ 



cos ax 



2 = 




cos [67 (4a)] \ (W\ f sin [o 2 /(4a)] 1 /b 2 ^ 
sin [o 2 /(4a)] J U ^4aJ \ cos [6 2 /(4a)] J ^4aJ 



(6) 



are differentiated under the integral sign we get the divergent integrals @. 
After some incorrect manipulations, Bierens de Haan obtains the value for 
these integrals |§, pp 443. He then differentiates under the integral sign to 
get for (|). 

The integrals in (H) and @ may be evaluated using the methods in the 
proof of the Proposition. 

The tables of Bierens de Haan had many errors, both mathematical and 
typographical. Two long works discussing the correctness of his tables are Eft 



and J^l . Neither mentions our divergent integrals. All of the integral tables 
listed above have received considerable scrutiny. The journal, Mathematics 
of computation, and its predecessor, Mathematical tables and other aids to 
computation, list numerous errata. However, despite over 300 published pages 
of errata related to the above tables there do not seem to be any references 
to ([!]) and (§). The article JT^| compares the correctness of various integral 



tables. The interested reader should consult this article to see how shockingly 
high the error rates are. 

4. Maple and Mathematica. Here are how Maple (V.4) and Mathemat- 
ica (4.0) fare. Maple correctly evaluates (||) and (|6|) for arbitrary a and b but 
falters when asked to perform the calculation with specific numerical values. 
For example, it gives J °° sin(3.1x 2 ) cos(2.2x) dx = 0. Maple correctly says 
the integrals (|l|) and (|2]) diverge. Mathematica fails in a different way. It 
correctly calculates (||) and (|]) (considerable simplification needed to obtain 
the form of @ above). But, Mathematica thinks ([1]) and (fj) converge! It 
gives the same incorrect values that are in the tables. 

5. Differentiation under the integral sign. Differentiating the conver- 
gent integrals (0) and (||) under the integral sign with respect to b yields the 
divergent integrals ([!]) and (Q). This doesn't mean the functions defined by 
d^) and (0) aren't differentiable it just means we cannot obtain their deriva- 
tives by differentiating under the integral. Could we have predicted this in 
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advance? This is a difficult problem. Suppose we have f(x,y) dy. A suf- 
ficient condition for differentiating under the integral, suitable for Riemann 
integrals, is that J fi(x,y) dy converge uniformly in x. See pEfl , pp 260. For 
Lebesgue integrals the dominating condition \ fi(x,y)\ < g(y) where g £ L 1 
suffices. For Riemann and Lebesgue integrals, necessary and sufficient con- 
ditions for differentiating under the integral sign are harder to come by. 
However, this is a much simpler problem when we use the Henstock integral. 
The solution depends on being able to integrate every derivative, a property 
not held by either the Riemann or the Lebesgue integral. The interested 



reader can see [[Jl|]. A good introduction to the Henstock integral is given in 

urn. 
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